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(57) Abstract 

The present invention relates to a method for estimating an impulse response and a receiver in a radio system where the signal to be 
sent comprises a known training sequence (208), which receiver comprises means (304) for sampling the received signal, and means (400) 
for calculating a first estimate for the impulse response by means of the known training sequence. To enable an accurate determination 
of impulse response, the receiver comprises means (406, 408) for making preliminary decisions of the received samples by means of the 
first impulse response estimate, means (410} for calculating an error value of the estimated samples and the received samples calculated by 
means of the preliminary decisions, means (410) for'calculating a second estimate of the impulse response by minimizing said error value, 
and means (410) for calculating a new estimaie for the impulse response by combining the first and second estimates to one another. 
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METHOD FOR ESTIMATING IMPULSE RESPONSE, AND RECEIVER 

FIELD OF THE INVENTION 

resoons* 7™ T*? lnVenti ° n re ' ateS t0 * f ° r SS,imatin 9 an im P^e 

response ,n a d.grtal radio system, where the signals to be transmitted 

5 compnse bursts formed of symbo«s, which bursts comprise a k loT^ g 

le'encl ' mPU ' Se by ^ ° f the fining 

BACKGROUND OF THE INVENTION 

10 station J« 3 ty ^ Ca ' r ,,U ' ar radi ° environfTient . *e signals between a base 

bit ' S Cnber termina ' eqU ' pment °n several routes 

between a transmitter and a receiver This mniti ™Z „ 
... lea(< hl< . , „ > e«-eiver. i ms multi-path propagation is mainly 

caused by Signal reflections from surrounding surfaces. Signa.* travelling on 

15 * . r ° U,eS a ' ** reCe " er at difte '«« *™ because of a d We re n" 

muTT y ' ™ S h ° ldS ^ fer •»* *•»"«» - Emission C 
mult-pa h propagation of a signal can be monitored in a receiver by measuring 

be n iZt e te respons : of me wce ^ si9nai ' -*» - ~ 

be,ng visible as peaks proportional to their signal strength. Figure 1 illustrates 

20 a 00 e a^T T UlSe rSSPOnSe * ™* « eMmpte ' ^ " » • ^rLoliTaS 
104 T^r* °' ^ " «" a «*- axis 102. Peaks 

104. 106 108 of the curve, indicate the strongest multi-path propagated 
components of the received signal. opagaiea 

a known , Pn ° r a " S °' UUOnS *" lm P ulse '«P°nse Is estimated by means of 

25 examT a 9 T"*™ ** 2 Urates by way of 

Z actu , °T GSM SyS,em """""^ « art a "« — b«3 200 

202, actual data ,n two parts 204, 206 and a. known training 203 sequence 

. ~ is%6 id b d : e °: \ buret in a «» - 

sequence ,s 26 bits. In known solutions, such as in the GSM system the 
30 rC e nr° nSe " eS " ma,ed ^ •» received ££Z 

L eZlr'"'" 9 T 6 " 08 ' 16 ^ * 3 26 * ,, - ton 9 " aW "3 «•*""- - used 

la ToZ T 2T response tep - When ^ ^ a,ity ° f ,h * ~— - 

methods ™ S w„' S * '° atta ' n 3 eStlmaU ° n ~'* **'*»«n 

methods. This will ,mp al r the performance of the receiver. 
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CHARACTERISTICS OF THE INVENTION 

The object of the present invention is to accomplish a more 
accurate method than the present methods for estimating an impulse 
response. 

5 This will be attained with a method as shown in the preamble, which 

is characterized by making preliminary decisions by means of the first impulse 
response estimate, calculating an error value of the estimated samples and 
the received samples calculated by means of the preliminary decisions, 
calculating a second estimate of the impulse response by minimizing said error 

10 value, calculating a new estimate for the impulse response by combining the 
first and second estimates to one another, and by using the estimate obtained 
in this way when estimating the received sequence again. The present 
invention also relates to a receiver in a digital radio system where the signals 
to be sent comprise bursts formed of signals, which bursts comprise a known 

15 training sequence, which receiver comprises means for sampling the received 
signal, and means for calculating a first estimate for the impulse response by 
means of the known training sequence. The receiver of the invention is 
characterized in that the receiver comprises means for making preliminary 
decisions of the received samples by means of the first impulse response 

20 estimate, means for calculating an error value of the estimated samples and 
the received samples calculated by means of the preliminary decisions, means 
for calculating a second estimate of the impulse response by minimizing said 
error value, means for calculating a new estimate for the impulse response by 
combining the first and second estimates to one another, and means for 

25 estimating the received sequence again by using the estimate obtained in this 
way. 

The method and "the receiver of the invention have several 
advantages in compiarison with prior art solutions. The accuracy of the impulse 
response estimate is greater in the solution of the invention than in prior art 

30 methods. Furthermore, if the speed of the terminal equipment is high, the 
training sequence in the middle of the burst cannot alone give a correct picture 
of the impulse response as the impulse response may change even during a 
burst. With the solution of the invention, it is possible to calculate an impulse 
response separately for the start portion of the burst and for the end portion of 

35 the burst. 



} 



WO 97/44916 PCT/FI97/00295 

3 

BRIEF DESCRIPTION OF THE FIGURES : 

In the following, the invention will be explained in more detail with 
reference to the examples of the appended drawings, wherein 

Figure 1 shows an example explained above of the impulse 
5 response of a received signal, 

Figure 2 shows a normal burst of the GSM system, 

Figure 3 illustrates the structure of the receiver of the invention by 
means of a block diagram, and 

Figure 4 illustrates in detail an example of the application of impulse 
10 response estimation. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The invention can be applied to any digital radio system where a 
burst comprises a training sequence. One example of this kind of a system is 
the GSM cellular radio system, and in the following, it is used as an example 

15 when explaining the invention without restricting thereto. 

In the solution of the invention, the impulse response of a channel is 
estimated adaptiveiy and in multi-phase. At first, a first preliminary estimate is 
calculated for the impulse response, by means of which, a preliminary bit 
decision is made. After this, an error value of the estimated samples and the 

20 , received samples calculated by means of the preliminary decisions is 
calculated, and a second estimation of the impulse response is calculated by 
minimizing said error value., In the preferred embodiment of the invention, the 
error value between the estimated samples and the received samples is 
calculated by forming first a convolution between the preliminary decisions and 

25 the first impulse response estimate and by summing the squared differences 
of the received samples and the estimated samples. 

After this, the first and second impulse response estimates are 
combined. In the preferred embodiment of the invention, the first and second 
impulse response estimates are .combined to one another by weighting the 

30 estimates by coefficients according to. their lengths and by, summing the 
received estimated tap values to one another. The actual detection will be 
carried out by using the impulse response obtained in this way. 

In the following, let us examine the structure of the receiver of the 
invention which is illustrated with relevant parts in the block diagram of Figure 

35 3. Both the base station and the subscriber terminal equipment may operate 
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- as receivers. of the invention. The receiver comprises an antenna 300, a signal 
received with which is conveyed to radio frequency parts 302 where the signal 
is converted to an intermediate frequency. From the radio parts, the signal is 
taken to converter means 304 where the signal is converted from an analog 
5 form into a digital form by sampling. A digital signal 306 propagates to means 
308 where the impulse response of the channel is estimated and detection is 
carried out by utilizing a Viterbi detector, for example. The detected signal is 
further conveyed to a channel decoder 312 and onwards 314 to other parts of 
the receiver. As is obvious to those skilled in the art, the receiver of the 

10 invention naturally comprises also other components than the ones described 
above, such as filters, but they are irrelevant to the invention and they have 
not been described for the sake of clarity. 

Let us ^examine next in more detail an example of the structure of 
the estimation block 308 which is illustrated with relevant parts in the block 

15 diagram of Figure 4. A received sampled signal 306 arrives at the estimation 
block as an input. The sampled signal is conveyed to a first estimation means 
400 where a first estimate of the impulse response is calculated with prior art 
methods by cross-correlating the samples with a training sequence. The first 
estimation block can be realized with methods known to those skilled in the 

20 art. 

A first estimate 416 of the impulse response is conveyed to a filter 
402 adapted by means of switches 412 and 414 and to a means 404 where 
autocorrelation taps are calculated. At this stage the switches 412 and 414 are 
in position 1. The received sampled signal 306 is also as an input of the 

25 adapted fitter. The adapted filter 402 restores the signal distorted in the 
channel into an original data stream by symbol error probability depending on 
interference factors, such as interference caused , by adjacent received brts. 
The autocorrelation taps of the estimated impulse response of the channel are 
formed *of impulse response information with a means 404.* The procedures 

30 described above can be realized by general and signal processors or by a 
detached logic circuit, for example. t , 

The outputs of the adapted filter 402 and the calculation means 404 
are conveyed to a prior art equalizer 406, preferably to a Viterbi equalizer. An 
output 310 of Viterbi 406 is conveyed to a decision means 408, where hard bit 

35 decisions are made on the basis of the first impulse response. 
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The output of the decision means 408 is coupled to an input of a 
second estimator 410. The received sampled signal 306 and the first estimate 
416 of the impulse response are also as an input of the second estimator 410. 
in the second estimator, a second estimate is calculated for the impulse 
5 response by means of the bit decisions and the received signal, and the first 
and second impulse response estimates are combined to one another. The 
second estimator can be realized by general and signal processors or by a 
detached logic circuit, for example. An estimate 418 obtained in this way is 
conveyed to the filter 402 adapted by means of the switches 412 and 414 and 

10 to the means 404 where the autocorrelation taps are calculated. At this stage 
the switches 412 and 414 are in position 2. On the basis of more exact 
impulse response information, the adapted filter 412 and the calculation 
means 414 calculate new values that are conveyed to the Viterbi equalizer 
406 from which the signal is taken to other parts of the receiver. 

15 The calculation of the second impulse response estimate of the 

invention is now examined in more detail. In the method of the invention, the 
squared error of samples is minimized by means of the received samples and 
the first impulse response estimate. The estimated samples are calculated as 
a convolution from hard bit decisions. It is assumed that a 5-bit impulse 

20 response is used, in which case the received bits are 

X k = (x kI x*_ 1f Xfc.j, x k _ 3 . x k _J T ; x € {-1 ,1} 
and the estimated impulse response is 

25 

W = (h 0 h, h 2 h 3 h 4 ) T . 
Therefore the estimated samples will b^ obtained now: 



30 Y = W 1 ^. r 

The error of the estimated samples will be determined next in 
comparison with the received samples. The' error can be calculated as a 
squared error, taking into account N bits at a time in the calculation, The error 
35 function can be formulated as follows: 
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* 2 = Z(-X*)-Y) 2 ,. 



where y(k) has received the samples and j is the beginning of a bit sequence 
being examined. A square of error, (MSE, minimum squared error) is thus 
used above as an error value. The error calculated in this way is minimized. 
5 The minimization can be carried out with various different methods described 
in the relevant literature, such as LMS t Kalman, direct matrix inversion, etc. 
For example, in the LMS method a gradient is taken from the function defined 
above and a solution is searched iteratively from the direction of the gradient. 
In principle, the square of the difference can be minimized by means of the 
10 mere difference without calculating the actual square. A new impulse response 
estimate will be obtained by means of the following formula: 



w 



e K = y(k)-Y 



where ji is a constant between (0 t 1). 

With value t=0. corresponding W 0 is the previously estimated first 
impulse response estimate W. This iterative formula can be repeated a desired 
number of times and after each repetition, an updated impulse response 
20 estimate can be obtained. It is assumed here that there are c times of 
repetitions, in which case the second impulse response estimate will thus be 
W c . 

In the following, the first and second impulse response estimates 
are combined by weighting the estimates by coefficients according to their 
25 lengths and by summing the corresponding estimated tap values with one 
another. Summing can be described with a formula 

m»W a + ;V»W c 

where m is the length of the bit sequence used in calculating the first impulse 
response estimate (that is, haff of a training sequence, 16 in the case of GSM) 
30 and N is the length of the bit sequence used in calculating the second impulse 
response estimate. 



# 
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The impulse response estimate W tot obtained in this way is used in 
the adapted fitter and in the calculation of autocorrelation taps, and final bit 
decisions are calculated by means of it ; 

It is also possible to calculate the impulse response estimate again 
5 a desired number of times by using the impulse response estimate W tot as an 
initial value when calculating new estimated samples and by minimizing the 
squared difference again. 

With the solution of the invention,, it is also possible to calculate an 
impulse response separately for the start portion of the burst and for the end 
10 portion of the burst. Then a more accurate result will be obtained especially if 
the speed of the terminal equipment is high. The first impulse response 
estimate is calculated by means of the training sequence in the middle of the 
burst, and the other impulse response estimates are calculated separately for 
the start portion of the burst and for the end portion of the burst 
15 Although the invention is explained above with reference to the 

example of the accompanying drawings, it is evident that the invention is not 
restricted thereto, but it can be modified in various ways within the inventive 
idea disclosed in the appended claims. 
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CLAIMS 

1. A method for estimating an impulse response in a digital radio 
system, where the signals to be transmitted comprise bursts formed of 
symbols, which bursts comprisie a known training sequence (208), in which 

5 method the received signal is sampled, and a first estimate (416) is calculated 
for the impulse response by means of the known training sequence, char- 
acterized by 

making preliminary decisions by means of the first impulse 
response estimate, 

10 calculating an error value of the estimated samples and the 

received samples calculated by means of the preliminary decisions, 

calculating a second estimate of the impulse response by 
minimizing said error value, 

calculating a new estimate (418) for the impulse response by 
15 combining the first and second estimates to one another, and by using the 
estimate obtained in this way when estimating the received sequence again. 

2. A method according to claim 1, characterized in that the 
second estimate is calculated rterath/ely a desired number of times again by 
using the new estimate obtained as an initial value the previous time when 

20 forming preliminary decisions. 

3. A method according to claim 1, characterized in that the 
first and second estimates are combined by weighting the estimates by 
coefficients according to their lengths and by summing the estimated tap 
values to one another. 

25 4. A method according to claim ^characterized in that the 

error value is calculated by forming first the estimated decisions by calculating 
a convolution between the preliminary decisions and the first impulse 
response estimate and by summing the squared differences of the received 
samples and the estimated samples. 

30 5. A method according to claim 1, characterized in that the 

second estimate is calculated separately for the start portion of the burst and 
for the end portion of the burst. 

6. A receiver in a digital radio system where the signals to be sent 
comprise bursts formed of signals, which bursts comprise a known training 

35 sequence (208), which receiver comprises means (304) for sampling the 
received signal, and means (400) for calculating a first estimate for the impulse 
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response by means of the known training sequence, characterized in 
that the receiver comprises 

means (406,408) for making preliminary decisions of the received 
samples by means of the first impulse response estimate, 
5 means (410) for calculating an error value of the estimated samples 

and the received samples calculated by means of the preliminary decisions, 

means (410) for calculating a second estimate of the impulse 
response by minimizing said error value, 

means (410) for calculating a new estimate for the impulse 
1 0 response by combining the first and second estimates to one another, 

and means (400) for estimating the received sequence again by 
using the estimate obtained in this way. 

7. A method according to claim 6, characterized in that the 
receiver comprises means (410) for combining the first and second estimates 

1 5 to one another by weighting the estimates by coefficients according to their 
lengths and by summing the corresponding estimated tap values to one 
another. 

8. A method according to claim 6, characterized in that the 
receiver comprises means (408, 410) for forming estimated decisions by~ 

20 calculating a convolution of the preliminary decisions and the firet impulse 
response estimate and by summing the squared difference of the received 
samples and the estimated samples. 
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